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I  INTRODUCTION 


The  relevance  of  naturally  occurring  F-region  plasma  instabilities 
for  understanding  the  onset,  the  nonlinear  evolution,  and  the  ultimate 
decay  of  nuclear-induced  strlatlons  is  well  established.  Indeed,  in 
lieu  of  actual  nuclear  data,  naturally  occurring  strlatlons  such  as 
those  in  the  auroral  and  equatorial  ionospheres  provide  the  only  means 
of  accumulating  large,  high-quality  data  bases  for  verifying  our  under¬ 
standing  of  these  processes;  moreover,  the  radio-wave  propagation 
effects  in  the  auroral  and  equatorial  regions  can  degrade  the  perfor¬ 
mance  of  DoO  systems. 

The  objective  of  the  work  conducted  under  this  contract  was  to  use 
our  extensive  auroral  zone  and  equatorial  data  base  to  identify  and  mea¬ 
sure  in  situ  the  parameters  that  affect  the  onset,  nonlinear  evolution, 
and  decay  of  strlatlons.  The  intensity  and  spectral  characteristics  of 
Intermediate-  and  small-scale  strlatlons  are  most  Important  because  they 
disrupt  radio-wave  systems. 

Investigations  of  equatorial  spread-F  (ESF)  phenomena  have  made 
valuable  contributions  to  this  effort.  Over  the  past  five  years,  the 
Defense  Nuclear  Agency  (DNA)  has  supported  a  number  of  field  campaigns 
that  have  used  extensively  the  ALTAIR  and  TRADEX  radars.  The  commonal¬ 
ity  of  the  physical  processes  that  structure  F-region  ionization  in  the 
equatorial  and  auroral  zones  (Ossakow,  1979)  has  been  thoroughly 
exploited.  The  equatorial  phenomena,  however,  are  much  more  consistent 
in  their  development  and  evolution,  evidently,  because  of  the  much 
broader  variety  (and  variability)  of  auroral-zone  phenomena  that  can 
affect  the  instability  development.  These  include  particle  precipita¬ 
tion,  a  highly  conducting  E  layer,  and  large  convection  velocities. 

In  Section  II  of  this  report,  we  review  the  high- latitude  work  that 
has  been  performed  under  this  contract.  This  work  includes  analysis  of 


5 


spaced-recelver  data  froto  the  Wideband  satellite  experiment  and  the 
development  of  a  theoretical  framework  for  Interpreting  the  global  mor¬ 
phology  of  high-latitude  scintillation.  The  model  has  been  used  to  for¬ 
mulate  a  set  of  hypotheses  that  can  be  tested  with  data  from  the  planned 
DNA  polar  satellite  experiment. 

Section  III  is  an  overview  of  all  the  experimental  ESF  work  that 
has  been  performed  under  DNA  programs.  This  contract  is  the  conclusion 
of  the  equatorial  experimental  program  for  hlgh-altltude  nuclear-effects 
predictions.  We  expect  that  future  equatorial  programs  will  use  the 
regular  occurrence  of  ESF  as  a  source  of  highly  structured  plasma  with 
known  spectral  characteristics  for  systems  tests,  or  testing  new  adapta¬ 
tions  of  the  basic  channel  modeling  method  for  evaluating  the  deleteri¬ 
ous  effects  of  propagation  disturbances,  or  both. 

In  Section  IV  we  discuss  the  ramifications  of  the  results  summa¬ 
rized  in  Sections  II  and  III  with  regard  to  probable  DoD  needs. 


II  HIGH-LATITUDE  STRUCTURE 


An  important  finding  from  the  Wideband  satellite  experiment  was  the 
nighttime  scintillation  enhancement  that  occurs  with  high  probability 
whenever  the  propagation  path  lies  within  the  plane  of  the  local  L  shell 
(Fremouw  et  al.,  1977;  Rlno  et  al».  1978).  The  phenomenon  has  been 
attributed  to  a  geometrical  enhancement  from  sheet-like  Irregularities. 
Evidence  Is  accumulating,  however,  for  a  localized  source  region  that 
also  contributes  to  the  nighttime  scintillation  enhancement.  Through 
combined  Incoherent  scatter  and  scintillation  observations,  the  F-region 
source  of  the  Irregularities  and  the  dynamics  of  the  source  region  have 
been  Identified  (Vickrey  et  al.,  1980).  Enhanced  scintillations  appear 
to  be  associated  with  the  horizontal  plasma -density  gradients  (or  walls) 
of  localized  plasma-density  enhancements,  or  (blobs),  in  the  auroral  F 
layer. 

The  experimental  and  theoretical  work  that  has  led  to  our  current 
understanding  of  this  hlgh-latltude  phenomenon  is  reviewed  in  Rlno  and 
Vickrey  (1982).  Because  of  the  exceptionally  good  definition  of  the 
source  region  and  Its  dynamics,  we  have  concentrated  our  analysis 
efforts  on  this  particular  source  of  hlgh-latltude  Irregularities. 

Under  active  conditions  In  the  nighttime  sector  of  the  auroral 
zone,  the  convection  electric  field  has  a  westward  component  (southward 
drift).  On  this  basis,  a  northward  directed  gradient  Is  stable  to  the 
gradient-drift  Instability  (e.g.,  Llnson  and  Workman,  1970).  This  led 
Ossakow  and  Chaturvedl  (1979)  to  propose  the  current-convective  insta¬ 
bility  as  the  structuring  mechanism.  Field-aligned  currents  provide  the 
free  energy  to  destabilize  the  gradients  that  are  otherwise  stable  to 
the  gradient-drift  Instability.  For  typical  observed  values  of  field- 
aligned  currents,  however,  this  Instability  Is  not  nearly  as  strong  as 
the  gradient-drift  Instability  (e.g.,  Kesklnen  and  Ossakow,  1982). 
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The  unstable  F-region  blob,  however,  often  has  a  much  more  complex 
structure  than  a  single  plasma  enhancement.  Figure  1  Illustrates  this 
and  shows  a  meridional  electron-density  contour  map  made  from  a 
Chatanlka  radar  elevation  scan.  Vickrey  et  al.  (1980)  have  shown  that 
scintillation  enhancements  occurred  when  the  features  marked  (1)  and  (2) 
In  the  figure  passed  through  the  propagation  path  of  the  TRIAD  satel¬ 
lite. 

The  main  blob  has  alternating  gradients  associated  with  a  large- 
scale  modulation  that  was  likely  to  have  been  imparted  by  the  soft- 
electron  precipitation  source  that  generated  the  plasma-density  enhance¬ 
ments.  From  this  figure  It  is  clear  that  some  segment  of  the  blob  would 
be  gradient-drift  unstable,  even  without  the  help  of  destabilizing 
field -aligned  currents. 

The  effects  of  the  gradient-drift  and  current-convective  instabili¬ 
ties  are  clearly  evident  In  a  recent  numerical  simulation  of  a  two-sided 
plasma  enhancement  (Kesklnen  and  Ossakow,  1982).  In  their  model  an 
electric  field  with  components  directed  northward  and  westward  are 
included,  together  with  a  typical  field-aligned  current.  Figure  2  shows 
the  plasma  structure  after  evolution  over  1600  s.  The  dominance  of 
gradient-drift  growth  over  current-convective  growth  is  clearly  evident. 
The  dominant  large-scale  structuring  occurred  on  the  northern  wall,  which 
Is  gradient-drift  unstable  because  of  the  presence  of  a  westward 
electric-field  component.  The  slight  tilt  in  the  interpenetrating  low- 
and  hlgh-denaity  regions  is  produced  by  velocity  shear  effects  induced 
by  the  presence  of  the  northward  electric  field  component  (e.g.,  Perkins 
and  Doles,  1975).  We  note  in  passing  that  the  large-scale  structure  is 
similar  to  those  produced  by  the  gravity-driven  Rayleigh-Taylor  insta¬ 
bility  in  the  equatorial  ionosphere  (see  Section  III). 

As  with  all  numerical  simulations  that  accurately  model  the  macro¬ 
scale  structures,  however,  the  above  type  of  simulation  does  not  retain 
sufficient  resolution  to  preserve  the  intermediate-  to  small-scale 
structures  that  cause  the  scintillation  enhancements.  It  is  important, 
nonetheless,  to  verify  that  the  structuring  process  is  initiated  at 
large  scales.  In  this  sense,  it  is  necessary  to  determine  whether  east- 
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FIGURE  1  MERIDIONAL  ELECTRON-DENSITY  PROFILE  SHOWING  STRUCTURED  F-REGION  PLASMA 
ENHANCEMENTS,  "BLOBS" 


EAST-WEST  DISTANCE  —  km 
Source:  Keskinan  and  Ossakow  [1982] 

FIGURE  2  SIMULATION  OF  THE  NONLINEAR  EVOLUTION  OF  CONVECTIVE 
INSTABILITY  FOR  AURORAL-ZONE  GEOMETRY 


west  structure  similar  to  those  in  Figure  2  are  indeed  observed  in  the 
auroral  ionosphere. 

Thus,  under  the  current  contract,  in  November  1981  a  series  of 
experiments  were  conducted  using  the  Chatanika  radar  to  map  the  large- 
scale  blob  structure  in  a  horizontal  plane  transverse  to  the  geomagnetic 
field.  This  was  accomplished  by  combining  the  data  from  several  rapid 
elevation  scans.  Figure  3  shows  the  preliminary  results  from  one  such 
measurement.  These  results,  obtained  in  the  local  midnight  sector, 
clearly  show  east-west  structure  along  the  southern  boundary  of  an  F- 
reglon  blob.  The  generally  northeast  orientation  of  the  mean  southern 
boundary  is  most  likely  caused  by  the  ionization  source.  The  large-scale 
irregularities  along  the  southern  boundary  was  consistent  with  gradient- 
drift  structuring  by  a  southward  blowing  neutral  wind.  (In  the  event  the 
east-west  electric  field  was  very  small.) 

Our  preliminary  results  show  that  we  have  good  semiquantitative 
agreement  with  the  predicted  nonlinear  evolution  of  auroral-zone  F- 
region-plasma  enhancements.  We  have  obtained  a  large  data  base,  and 
detailed  analysis  is  being  pursued  under  a  separate  contract. 

Turning  now  to  the  intermediate-  to  small-scale  structure,  we  have 
performed  a  detailed  analysis  of  the  irregularity  anisotropy  as  measured 
by  the  Wideband  satellite  spaced-receiver  system  operated  at  Poker  Flat, 
Alaska.  The  measurement  technique  is  described  in  Rlno  and  Livingston 
(1982);  the  results  of  the  analysis  are  described  in  Livingston  et  al. 
(1982). 

By  carefully  measuring  the  variations  of  the  axial  ratio  and  orien¬ 
tation  angle  of  the  diffraction  pattern  during  a  large  number  of  satel¬ 
lite  passes,  Livingston  et  al.  (1982)  were  able  to  identify  systematic 
transitions  among  three  generic  types  of  irregularity  structures:  rods, 
wings,  and  sheets.  Rods  are  isotropic  in  the  plane  transverse  to  the 
geomagnetic  field  and  can  be  characterized  by  an  axial  ratio  (parallel 
to  transverse  dimensions)  of  a:l.  Wings  and  sheets  are  elongated  both 
along  the  geomagnetic  field  and  in  the  transverse  plane  along  the  direc¬ 
tion  of  the  local  L  shell;  hence,  they  are  characterized  by  two  axial 
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ratio  parameters  as  a:b:l.  For  sheets,  a  *  b.  Wings  are  intermediate- 
scale  structures  with  b  *  a/ 2. 

Figure  4  shows  the  general  pattern  established  by  Livingston  et  al. 
(1982).  Sheets  and  wings  are  confined  to  the  southern  region  of  the 
scintillation  zone.  The  axial  ratio  of  the  rods  systematically 
decreases  toward  the  poleward  boundary  of  the  scintillation  zone.  Wings 
are  generally  confined  to  the  near  midnight  sector.  This  clearly  sug¬ 
gests  a  close  connection  with  the  general  pattern  of  nighttime  auroral- 
zone  convection,  as  shown  schematically  in  Figure  5.  These  results  and 
their  ramifications  are  discussed  in  detail  in  Livingston  et  al.  (1982). 

To  provide  a  theoretical  framework  for  interpreting  the  morphology 
of  scintillation-producing  irregularities,  the  effects  of  cross-field 
diffusion  in  the  presence  of  a  conducting  E  layer  have  been  analyzed 
(Vickrey  and  Kelley,  1982).  By  postulating  a  source  region  at  the 
equatorward  boundary  of  the  polar  cap  and  a  detailed  convection  model, 
the  lifetime  of  kilometer-scale  Irregularities  was  computed.  The 
details  of  the  model  and  its  ramifications  are  discussed  by  Vickrey 
(1982). 

The  general  features  of  the  model  are  Illustrated  in  Figure  6, 
which  shows  the  diffusive  decay  of  a  kilometer-scale  structure  as  it 
convects  through  the  polar  cap  and  into  the  auroral  zone.  As  discussed 
by  Vickrey  (1981),  the  model  reproduces  the  general  features  of  morpho¬ 
logical  studies  of  hlgh-latltude  irregularities. 

To  summarize,  we  have  measured  the  east-west  structure  that  is  pre¬ 
dicted  theoretically  by  numerical  simulations  of  the  nonlinear  evolution 
of  a  two-sided,  F-region  plasma  enhancement  in  the  nighttime  auroral 
zone.  We  have  also  established  a  distinct  latitudinal  and  magnetic- 
local-time  variation  in  the  anisotropy  of  nighttime  auroral-zone  irregu¬ 
larities.  Finally,  we  are  developing  a  model  that  contains,  at  least  in 
rudimentary  form,  all  the  essential  elements  that  must  be  included  in  a 
truly  predictive  (rather  than  average)  hlgh-latltude  scintillation 
model. 


FIGURE  5  SCHEMATIC  DIAGRAM  OF  ANISOTROPY  VARIATION  AND  ITS  RELATION  TO 
AVERAGE  NIGHTTIME  CONVECTION  PATTERN 
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FIGURE  6 


MIDNIGHT 

MODEL  COMPUTATION  OF  DIFFUSIVE  DECAY  OF  KILOMETER 
SCALE  IRREGULARITIES 


The  same  considerations  apply  to  the  nonlinear  evolution  of  stria- 
tlons  that  are  generated  by  high-altitude  nuclear  detonations.  In  an 
early-time  nuclear  environment,  of  course,  conditions  are  radically  dif 
ferent  from  anything  that  occurs  in  the  natural  ionosphere.  The  very 
long  lifetimes  of  F-reglon  plasma  structures  suggest,  however,  that  the 
basic  features  of  the  model  being  developed  at  SRI  will  ultimately  domi 
nate  the  distribution  and  configuration  of  structure.  There  is  much 
work  to  be  done,  and  the  planned  DNA  High  Latitude  satellite  will  pro- 
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Ill  EQUATORIAL  STRUCTURE 


In  1977  DNA  initiated  an  investigation  of  plasma-density  structure 
in  the  nighttime  equatorial  ionosphere.  Interest  in  equatorial  struc¬ 
ture,  collectively  called  ESF,  was  instigated  on  the  hypothesis  that  the 
collisional  Rayleigh-Taylor  instability  was  the  dominant  structuring 
mechanism.  This  instability  is  the  gravitational  analog  of  the 
gradient-drift  instability,  the  principal  structuring  mechanism  in  the 
late-time  nuclear  environment.  Investigations  of  ESF  phenomena  were 
conducted  from  the  Kwajalein  Atoll,  Marshall  Islands,  in  field  experi¬ 
ments  during  each  of  five  summers,  1977  through  1981,  including  a  highly 
successful  rocket  (Plumex  I)  campaign  in*^'1979. 

In  the  following  subsections,  we  summarize  the  experimental  results 
obtained  from  this  series  of  field  programs  and  discuss  what  we  consider 
to  be  the  remaining  major  questions.  The  results  are  divided  into  four 
categories,  according  to  irregularity  scale  size:  large,  intermediate, 
transitional,  and  small. 


A.  Large-Scale  Structure  (Greater  than  10  km) 

Investigations  of  ESF  irregularities  have  shown  that  the  structur¬ 
ing  process  is  Initiated  at  large  scales  (l.e.,  greater  than  10  km)  and 
proceeds  through  an  apparent  hierarchy  of  plasma  Instability  processes 
that  produce  smaller  scales  (e.g.,  Haerendel ,  1973).  The  structuring 
mechanism  acting  at  these  largest  scales  appears  to  be  the  collisional 
Rayleigh-Taylor  instability.  During  the  early  stages  of  theoretical 
development,  the  collisional  Rayleigh-Taylor  instability  was  proposed  as 
a  possible  mechanism  (Dungey ,  1956).  However,  the  theory  was  played 
down  because  irregularity  growth  appeared  to  be  possible  only  in  the 
bottomside  F  layer,  whereas  ESF  also  occurred  in  the  topside  F  layer 
(e.g.,  Farley  et  al.,  1970).  Scannapieco  and  Ossakow  (1976)  have  since 
shown  that  the  nonlinear  collisional  Rayleigh-Taylor  instability  could 


account  for  irregularities  found  both  in  the  bottomside  and  topside  of 
the  F  layer. 

The  Kwajalein  experimental  results  have  shown  conclusively  that  the 
nonlinear  numerical  simulations  of  the  collislonal  Rayleigh-Taylor 
instability  (Scannapieco  and  Ossakow,  1976;  Ossakow  et  al.,  1979; 

Zalesak  and  Ossakow,  1980)  correctly  describe  the  development  of  large- 
scale  plasma  depletions,  or  bubbles.  The  procedure  for  experimental 
verification  was  as  follows.  First,  radar  backscatter  plumes,  observed 
with  the  ALTAIR  radar  (Tsunoda  et  al.,  1979)  were  shown  to  be  collocated 
with  plasma  bubbles  (Tsunoda  and  Towle,  1979;  Tsunoda,  1980a, b;  Towle , 
1980;  Szuszczewicz  et  al.,  1980).  The  temporal  evolution  of  plasma  bub¬ 
bles  could  then  be  followed  by  examining  a  time  sequence  of  spatial  maps 
of  ESF  backscatter  obtained  with  ALTAIR,  a  steerable  radar. 

Primary  structuring  of  large-scale  irregularities  began  in  the 
bottomside  F  layer  in  the  form  of  upwelllngs  (Tsunoda,  1981;  Tsunoda  and 
White,  1981).  Examples  of  bottomside  upwelllngs  can  be  seen  in  Figure 
7.  Three  upwelllngs  appear  in  the  figure  as  altitude  modulations  in  the 
isodensity  contours.  The  upwelllngs  then  push  upward,  forming  plasma 
bubbles,  and  penetrate  into  the  topside  F  layer.  The  upward  development 
of  plasma  bubbles,  mapped  spatially  as  radar  backscatter  plumes,  is 
Illustrated  in  Figure  8.  The  growth  of  plasma  bubbles  into  vertically 
elongated  "wedges"  was  similar  to  the  numerical  simulations  of  the  col- 
lisional  Rayleigh-Taylor  instability  by  comparing  the  spatial  relation¬ 
ship  of  backscatter  plumes  to  plasma  depletions  measured  in  situ  in  the 
neck  region  of  plumes  by  the  Atmospheric  Explorer  E  satellite  (Tsunoda 
et  al..  1982). 

The  description  of  large-scale  structure  (i.e.,  plasma  bubbles)  as 
a  two-dimensional  process  appears  justified  provided  that  the  descrip¬ 
tion  uses  flux-tube-integrated  quantities  (Haerendel,  1973;  Anderson  and 
Haerendel ,  1979).  The  field-aligned  nature  of  plasma  bubbles  was  sub¬ 
stantiated  with  ALTAIR  by  mapping  a  plasma  bubble  (using  incoherent- 
scatter  measurements)  over  more  than  ten  degrees  of  magnetic  dip  latitude 
(Tsunoda,  1980b).  A  map  of  a  plasma  bubble  in  altitude  and  magnetic-dip 
latitude  is  shown  in  Figure  9. 
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BEFORE  THE  DEVELOPMENT  OF  ESF  BACKSCATTER 


In  addition  to  the  primary  structuring,  ALTAIR  results  Indicated 
that  secondary  structuring  occurred  regularly  along  the  western  walls  of 
bottomslde  F-layer  upwelllngs  (Tsunoda.  1981;  Tsunoda  and  White,  1981). 


Tsunoda  (1981)  suggested  that  this  secondary  structuring  was  produced  by 
the  gradient-drift  instability  driven  ly  an  eastward  neutral  wind  blow¬ 
ing  through  the  western  walls  of  upwelllngs.  In  Figure  10,  the 
structuring  of  the  western  wall  of  an  upwelling  is  illustrated.  By 
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tracking  the  plume  seen  in  the  upper  left  hand  corner  of  Figure  8,  we 
find  that  the  western  wall  continues  to  structure  over  a  period  of  more 
than  an  hour,  even  after  the  decay  of  the  primary  plume.  The  structur¬ 
ing  appears  in  the  form  of  smaller-scale  plasma  bubbles  or  irregulari¬ 
ties.  The  strong  backscatter  from  the  western  wall  and  the  absence  of 
backscatter  from  the  eastern  wall  strongly  supports  the  proposed 
gradient-drift  mechanism.  (We  should  note,  of  course,  that  because  the 
upwelling  is  located  in  the  bottomside  F  layer,  the  collisional 
Rayleigh-Taylor  instability  is  also  operative  there.)  Szuszczewicz  et 
al.  (1981)  showed  that  the  Plumex  II  rocket  measurements  (launched 
0957:30  UT)  found  relatively  weak  plasma -density  structure  within  the 
eastern  wall  of  the  upwelling  in  Figure  10(e). 

Conclusive  evidence  for  the  action  of  the  gradient-drift  instabil¬ 
ity  along  the  western  wall  of  the  upwelling  is  presented  in  Figure  11. 
The  east-west  plasma  drift  was  estimated  by  measuring  the  horizontal 
displacements  of  traceable  backscatter  features  in  Figure  8.  Tsunoda  et 
al.  (1981)  showed  that  the  east-west  plasma  drift  is  characterized  by  a 
velocity  shear  with  altitude.  The  velocity  shear  is  evident  in  Figure 
11  by  noting  that  the  solid  bars  correspond  to  high-altitude  measure¬ 
ments  and  the  dashed  bars  correspond  to  low-altitude  measurements.  The 
velocity  difference  exceeds  100  m/s.  Figure  11  also  shows  the  eastward 
neutral  wind,  measured  simultaneously  with  a  Fabry-Perot  interferometer 
(Sipler  et  al. ,  1980).  Fabry-Perot  interferometer  measurements  refer  to 
the  bottomside  of  the  F  layer  in  which  630.0-nm  emissions  are  most 
Intense.  The  neutral  wind  is  seen  to  accelerate,  starting  around  0745 
UT,  and  approach  the  mean  eastward  plasma  drift  at  high  altitudes.  The 
ion-neutral  slip  velocity  that  drives  the  gradient-drift  instability  is 
the  difference  between  the  low-altitude  plasma  drift  (dashed  bars)  and 
the  neutral  wind.  From  this  figure,  there  is  little  doubt  that  the 
western  wall  was  structured  primarily  by  the  gradient-drift  instability. 
Figure  12  summarizes  the  experimental  results  on  large-scale  ESF  struc¬ 
ture,  published  to  date.  Structuring  appears  to  be  initiated  by  an  up¬ 
welling  in  the  bottomside  F  layer  [Figure  12(a)].  The  upwelling  pushes 
upwards  [Figure  12(b)]  via  the  collisional  Rayleigh-Taylor  instability. 


(e) 
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FIGURE  12  DESCRIPTIVE  MODEL  OF  PLASMA  STRUCTURING  IN  THE  EQUATORIAL 
F  LAYER  PRODUCED  BY  THE  GRADIENT-DRIFT  AND  COLLISIONAL 
RAYLEIGH-TAYLOR  INSTABILITIES 
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producing  a  plasma  bubble  that  penetrates  into  the  topside  F  layer 
[Figure  12(c)].  The  bubble  continues  to  develop  upward  until  it  reaches 
an  altitude  at  which  the  ambient  plasma  density  is  identical  to  that 
within  the  bubble  [Figure  12(d)].  Structuring  of  the  eastern  wall  of 
the  wedge  (formed  by  the  upward-developing  plasma  bubble)  at  high  alti¬ 
tudes  might  occur  via  the  Rayleigh-Taylor  instability,  which,  together 
with  cross-field  diffusion,  could  Chen  act  to  close  Che  wedge  [Figure 
12(e)].  Ac  Che  same  time,  Che  western  wall  of  the  upwelling  (in  Che 
bottomslde  F  layer)  begins  to  structure  via  the  gradient-drift  instabil¬ 
ity,  driven  by  an  eastward  neutral  wind  [Figures  12(d)  and  12(e)].  The 
slip  velocity  is  enhanced  in  the  upwelling  by  a  velocity  shear  that  is 
produced  by  the  F-region  dynamo  acting  in  the  presence  of  an  altitude- 
dependent  Pedersen  conductivity.  The  analogy  of  the  western  wall  of  the 
upwelling  to  the  well-known  structuring  of  a  barium  ion  cloud  is  illus¬ 
trated  in  Figure  12(f). 

More  recent  (and  as  yet  unpublished)  results  of  ALTAIR  data  analy¬ 
sis  indicate  Chat  the  scenario  depicted  in  Figure  12  is  not  the  only 
process  by  which  plasma  bubbles  and  backscatter  plumes  develop.  Figure 
13  illustrates  another  scenario.  The  sequence  of  contour  maps  con¬ 
structed  from  ALTAIR  backscatter  data  taken  on  22  July  1979  reveals  bub¬ 
ble  and  plume  development  without  the  presence  of  a  large-scale  upwell¬ 
ing.  In  Figure  13(a),  we  see  Che  occurrence  of  a  wedge  or  bubble  in  the 
bottomslde  F  layer  with  backscatter  associated  with  its  upper  wall. 

There  is  no  evidence  of  any  large-scale  wave  structure  in  the  bottomslde 
plasma-density  contour,  as  reported  by  Tsunoda  and  White  (1981).  [For 
example,  compare  Figure  13(a)  to  Figure  7  or  8(a)]. 

In  Che  next  three  maps  [Figures  13(c)  to  13(e)],  we  observe  a 
slight  clockwise  rotation  of  the  plume  as  it  continues  to  develop  up¬ 
ward.  Moreover,  a  "kink"  develops  in  the  neck  of  the  plume  [Figure 
13(d)].  Both  the  kink  and  the  clockwise  rotation  can  be  explained  in 
terms  of  a  velocity  shear.  The  plasma-density  profile  taken  along  the 
radial  line  to  the  right  of  the  plume  in  Figure  13(e),  is  shown  in 
Figure  14.  We  find  that  the  peak  of  the  F  layer  occurred  around  the 
S30-kffl  altitude  and  the  steep  bottomslde  gradient  occurred  below  about 
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FIGURE  13  SEQUENCE  OF  ALTAIR  BACKSCATTER  MAPS 
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FIGURE  14  ELECTRON-DENSITY  PROFILE  OBTAINED 
ALONG  RADIAL  LINE  IN  FIGURE  13(e) 


400  km.  The  clockwise  rotation  of  the  plume,  therefore,  appears  to  be 
associated  with  the  slowly  varying  portion  of  the  plasma -density  profile 
and  the  kink  with  the  steep  bottomside  profile.  This  interpretation  is 
consistent  with  the  numerical  simulations  of  Rayleigh-Taylor  bubbles  by 
Zalesak  et  al.  (1982).  In  their  simulation,  they  showed  that  in  the 
presence  of  finitely  conducting  E  layer,  the  velocity  shear  is  directly 
related  to  the  F-layer  Pedersen  conductivity  profile. 

A  similar  plume  development  sequence  is  seen  in  the  contour  maps 
presented  in  Figure  15.  In  these  maps,  which  follow  immediately  after 
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the  last  map  in  Figure  13,  we  see  that  appearance  of  another  wedge 
[Figure  15(a) J,  followed  by  its  development  into  a  small  plume  in 
Figures  15(c)  to  15(e).  Similar  velocity  shear  effects  are  also 
present.  It  is  interesting  to  note  that  in  spite  of  the  initial  absence 
of  large-scale  wave  structure  and  upwellings  in  the  bottomside  F  layer, 
the  plumes  are  still  spaced  a  few  hundred  kilometers  apart;  in  this 
case,  350  km.  There  is  also  evidence  in  Figure  15  that  a  large-scale 
wave  structure  (upwellings)  appears  to  develop  in  Figures  15(c)  to 
15(e),  after  the  development  of  at  least  two  of  the  plumes.  The  plume 
near  the  left  side  of  all  the  maps  in  Figure  15,  however,  might  have 
been  generated  after  the  development  of  upwellings. 

The  above-described  results  indicate  that  large-scale  wave  struc¬ 
ture  does  not  have  to  be  present  to  initiate  the  generation  of  bubbles 
and  plumes.  Yet  Tsunoda  and  White  (1981)  showed  that  the  collisional 
Rayleigh-Taylor  instability  does  not  appear  to  be  capable  of  amplifying 
thermal  fluctuations  to  detectable  levels,  and  that  a  seed  perturbation 
appears  necessary. 

The  question  as  to  what  controls  bubble  or  plume  development  remains. 
This  question  is  analogous  to  what  controls  the  bifurcation  tendency  of 
barium  (or  nuclear-produced)  plasma  clouds.  In  the  case  of  the  initial 
formation  of  bubbles,  we  are  dealing  with  a  slab-like  geometry,  although 
upwellings  provide  a  more  concave  surface.  In  the  case  of  bifurcation  of 
plasma  bubbles  or  plumes,  we  are  dealing  with  a  much  smaller  radius  of 
curvature. 

An  example  of  plume  (and  bubble)  bifurcation  can  be  seen  in  the 
plume  in  Figure  10(a).  There  is  a  distinct  indentation  in  the  head  of 
the  plume,  apparently  forming  two  separate  secondary  bubbles.  All 
plumes,  however,  do  not  bifurcate  to  the  same  extent.  An  example  is  the 
left  plume  in  Figure  15(e).  A  later  map  of  the  same  plume  is  presented 
in  Figure  16.  There  is  no  obvious  bifurcation.  Finally,  we  present  an 
outstanding  example  of  extensive  plume  bifurcation  in  Figure  17.  This 
plume,  observed  on  10  July  1980,  displays  a  two-step  bifurcation  his¬ 
tory.  Apparently,  the  primary  bubble  (plume)  first  trifurcated  to  form 
secondary  bubbles  labeled  1,  2,  and  3.  Secondary  bubble  1  then  trifur¬ 
cated  to  form  tertiary  bubbles  labeled  a,  b,  and  c. 
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FIGURE  15  CONTINUATION  OF  THE  SEQUENCE  OF 
ALTAI R  BACKSCATTER  MAPS  SHOWN  IN 
FIGURE  13 
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FIGURE  17  ALTAIR  BACKSCATTER  MAP  SHOWING 
PROGRESSIVE  STRUCTURING  OF  BACK 
SCATTER  PLUME 


B.  Intermediate-Scale  Structure  (10  km  to  100  m) 

The  spatial  wavelengths  labeled  as  Intermediate  scale  structure 
include  the  short  end  of  the  irregularity  spectrum  believed  to  be  pro¬ 
duced  by  the  colllslonal  Rayleigh-Taylor  instability.  At  a  100-m  spa¬ 
tial  wavelength,  finite  gyroradlus  effects  are  still  minimal,  particu¬ 
larly  in  a  colllslonal  plasma  (Hudson  and  Kennel,  1975).  At  shorter 
wavelengths  (transitional  and  small-scale  structure),  other  plasma  pro¬ 
cesses  begin  to  dominate  the  production  of  structure. 

Our  understanding  of  intermediate-scale  structure  has  been  advanced 
by  scintillation  measurements  and  by  in  situ  satellite  and  rocket  mea¬ 
surements.  Scintillation  measurements  were  made  at  Kwajalein  by  using 
the  polar-orbiting  Wideband  satellite  and  the  Plumex  I  rocket.  In  situ 
measurements  also  were  made  by  the  Plumex  I  rocket. 

Phase-scintillation  measurements  using  the  Plumex  I  rocket  beacon 
contain  evidence  that  the  principal  scintillation-producing  irregulari¬ 
ties  are  Imbedded  in  high  plasma-density  regions  adjacent  to  plasma  bub¬ 
bles  (Rino  et  al. ,  1981).  Strong  scintillations  were  not  associated 
with  plasma-depleted  regions.  This  conclusion,  however,  requires 
further  verification  because  the  plasma  bubble  probed  by  the  Plumex  1 
rocket  was  decaying.  It  is  possible  that  intermediate-scale  structure 
is  associated  also  with  plasma-depleted  regions  during  the  growth  phase 
of  plasma  bubbles. 

Perhaps  a  more  significant  finding,  from  in  situ  rocket  and  scin¬ 
tillation  measurements,  was  that  the  spectral  slope  associated  with 
intermediate-scale  structure  depends  on  the  strength  of  the  irregulari¬ 
ties  (Livingston  et  al.,  1981;  Rino  et  al.,  1981).  Both  Wideband  and 
rocket  beacon-scintillation  measurements  and  in  situ  AE-E  satellite  mea¬ 
surements  contain  evidence  that  the  slope  tended  to  become  more  shallow 
with  increasing  turbulence  strength.  As  yet,  there  is  no  satisfactory 
explanation  for  the  slope  dependence  on  turbulence  strength. 
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In  addition  to  depending  on  turbulence  strength,  the  spatial  irreg¬ 
ularity  spectrum  measured  by  the  Plumex  I  rocket  was  best  represented  by 
two  spectral  slopes,  with  the  breakpoint  located  around  500  m  (Rino  et 
al. ,  1981).  An  example  of  the  irregularity  spectrum  obtained  during  the 
rocket  up leg  at  the  316-km  altitude  is  shown  in  Figure  18.  This  dual 
slope  characteristic  appeared  confined  to  lower  altitudes  at  which  the 
largest  perturbations  developed  as  shown  in  Figure  19.  The  spectral 
slopes  have  been  plotted  as  a  function  of  altitude.  The  two  spectral 
Indices  converge  towards  a  value  near  two  at  an  altitude  between  350  and 
400  km.  Rino  et  al.  (1981)  suggested  that  the  spatial  wavelength  of 
500  m  might  be  a  manifestation  of  the  stirring  scale  length  associated 
with  two-dimensional  plasma  turbulence,  e.g.,  in  the  scenario  proposed 
by  Kelley  and  Ott  (1978). 

An  alternate  explanation  for  variations  in  spectral  slope  would  be 
the  presence  of  more  than  a  single  driver.  If  the  outer  scale  of  the 
irregularity  spectrum  is  determined  by  the  initial  gradient  scale  length 
(Kesklnen  et  al.,  1980),  two  drivers  might  be  acting  on  different  ini¬ 
tial  gradient  scale  lengths  to  produce  a  composite  spectrum  that  is 
characterized  by  two  spectral  slopes.  Tsunoda  (1981)  proposed  that,  in 
addition  to  the  colllslonal  Rayleigh-Taylor  instability,  the  gradient- 
drift  instability  was  probably  operative.  With  the  gradient-drift 
instability  acting  on  the  western  wall  of  large-scale  upwellings  (see 
Section  III-A),  driven  by  an  eastward  neutral  wind,  we  might  expect  an 
input  gradient  scale  length  associated  with  the  western  wall  or  with 
primary  irregularities.  Because  the  velocity  shear  that  allows  the 
gradient-drift  instability  to  be  operative  is  present  only  in  the 
bottomslde  F  layer,  it  is  tempting  to  associate  the  altitude  dependence 
seen  in  Figure  19  with  a  dual  driver  mechanism.  This  concept,  however, 
still  remains  to  be  clarified. 

C.  Transitional-Scale  Structure  (100  m  to  10  m) 


The  transitional-scale  regime  represents  that  portion  of  the  spa¬ 
tial  irregularity  spectrum  that  is  produced  by  decreasing  contributions 
by  the  collisional  Rayleigh-Taylor  instability,  and  increasing  contribu- 
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FIGURE  19  MEASUREMENTS  OF  POWER-LAW  INDEX  IN  LARGE-SCALE  AND 
SMALL-SCALE  SPECTRAL  REGIMES  FROM  DNA  PLUMEX 
ROCKET  CAMPAIGN 


tions  by  other  plasma  processes.  The  collisional  Rayleigh-Taylor  insta¬ 
bility  is  damped  at  these  wavelengths  by  finite  gyroradius  effects. 

Other  plasma  processes  include  excitation  of  drift  waves  and  diffusion. 
Drift-wave  instabilities  are  presumably  driven  by  the  gradients  associ¬ 
ated  with  larger-scale  irregularities  produced  by  the  collisional 
Rayleigh-Taylor  instability.  It  is  not  yet  clear  whether  linear  excita¬ 
tion  of  drift  waves  occurs  or  whether  turbulent  cascade  occurs  through 
wave-wave  coupling. 

Structure  at  these  wavelengths  was  measured  in  situ  by  plasma 
probes  on  the  Plumex  I  rocket.  These  structures  correlated  with  the 
intermediate-scale  structure.  This  correlation  is  consistent  with  the 
concept  that  low-frequency  drift  waves  are  driven  by  the  larger-scale 
irregularities. 

Kelley  et  al.  (1982)  showed  that  the  low-altitude  spectrum  differed 
from  the  hlgh-altltude  spectrum  in  that  the  former  could  be  expressed 
with  two  spectral  slopes  whereas  the  latter  was  best  described  with  a 
single  spectral  slope.  The  spectral  slope  varied  from  -2  to  -5  with  a 
changeover  occurring  around  an  irregularity  wavelength  of  22  m.  The 
high-altitude  spectrum  was  characterized  by  a  spectral  slope  of  -5. 
Kelley  et  al.  (1982)  thought  that  the  -5  slope  was  consistent  with 
drift  waves  and  that  the  crossover  wavelength  was  the  ion  gyroradius. 

The  confinement  of  the  two-slope  spectrum  to  low  altitudes  was  inter¬ 
preted  in  terms  of  an  altitude  threshold  for  drift-wave  turbulence. 

That  is,  anomalous  diffusion  was  assumed  operative  above  the  280-km 
altitude,  and  collisional  diffusion  was  assumed  operative  below  that 
altitude. 

D.  Small-Scale  Structure  (Less  than  10  m) 


Small-scale  structure  includes  all  irregularities  with  spatial 
wavelengths  less  than  the  Ion  gyroradius.  In  this  regime,  the  colli¬ 
sional  Rayleigh-Taylor  (not  the  gradient-drif t)  Instability  is  applica¬ 
ble.  Even  the  low-frequency  collisional  drift  waves  are  strongly  damped 
by  ion  viscosity  (Huba  and  Ossakow,  1979).  This  regime  is,  therefore, 
believed  a  product  of  a  turbulent  cascade  process  in  which  wave-wave 


coupling  is  important,  and  on  which  the  high-frequency  drift  wave  insta¬ 
bility  might  be  operative. 

An  Important  aspect  of  the  Kwajalein  equatorial  program  was  the 
investigation  of  small-scale  irregularities  using  backscatter  radars. 
Backscatter  measurements  of  1-m  and  36-cm  irregularities  with  the  ALTAIR 
radar  (e.g.,  Tsunoda  et  al.,  1979;  Towle ,  1980)  were  complemented  by  11- 
cm  irregularity  measurements  with  the  TRADEX  radar  (Tsunoda ,  1980c). 
These  radar  investigations  provided  characterization  of  both  the  large- 
scale  structural  evolution  in  space  and  time  (Section  III-A)  and  small- 
scale  irregularity  growth  and  decay. 

Initial  ALTAIR  measurements  of  ESP  backscatter  in  1977  led  to  the 
realization  that  these  ESP  Irregularities,  with  scale  sizes  smaller  than 
the  ion  gyroradius,  could  not  be  produced  by  low-frequency  drift-wave 
instabilities  [Hudson  and  Kennel,  1975;  Costa  and  Kelley,  1978(a,b)J. 
Huba  et  al.  (1978)  proposed  that  the  observed  small-scale  structure  is 
likely  produced  by  a  high-frequency  drift-wave  instability,  in  particu¬ 
lar,  the  lower-hybrid-drift  (LHD)  Instability.  The  plasma-density  gra¬ 
dients  are  presumed  present  at  the  large-scale  end  of  the  irregularity 
spectrum  by  the  collisional  Rayleigh-Taylor  instability. 

The  LHD  instability  appeared  to  be  linearly  unstable  under  presum¬ 
ably  typical  ESP  conditions.  In  a  sample  calculation,  Huba  et  al. 

(1978)  estimated  that  a  gradient  scale  length,  L,  of  73  m  was  required 
for  instability  growth.  Costa  and  Kelley  (1978)  have  reported  in  situ 
measurements  of  L  less  than  75  m.  Moreover,  the  large  growth  rates  of 
this  instability  are  consistent  with  observations  of  sudden  onsets  of 
ESP  backscatter. 

Stabilizing  mechanisms  for  the  LHD  instability  have  since  been  con¬ 
sidered  and  found  to  impose  a  plasma-density  and  altitude  threshold  on 
small-scale  ESP  irregularity  growth.  Sperling  and  Goldman  (1980) 
pointed  out  the  Importance  of  electron-ion  collisions  as  a  damping  mech¬ 
anism.  To  Illustrate  this  effect,  Sperling  and  Goldman  (1980)  used  a 
plasma  density  of  10^  el/ca?  and  the  plasma  description  used  by  Huba  et 
al.  (1978)  and  showed  that  the  LHD  instability  was  linearly  stable. 
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They,  of  course,  pointed  out  that  instability  was  possible  in  low- 
plasma-density  regions,  such  as  in  plasma  bubbles.  Huba  and  Ossakow 
(1981)  also  noted  that  electron-neutral  collisions  will  damp  irregular¬ 
ity  growth  at  lower  altitudes  (below  240  km). 

The  general  preference  of  ESF  backscatter  for  low-plasma-density 
regions  has  been  experimentally  verified.  The  spatial  coincidence  of 
ESF  backscatter  with  low-plasma-denslty  regions  has  been  shown  through 
comparisons  with  total  electron  content  measurements  (Tsunoda  and  Towle, 
1979;  Rino  et  al . ,  1981),  incoherent-scatter  radar  measurements 
(Tsunoda,  1980a, b,  1981;  Towle ,  1980),  and  in  situ  measurements 
(Szuszczewicz  et  al.,  1980;  Kelley  et  al..  1982;  Tsunoda  et  al.,  1982). 

A  preference  also  for  upward-directed,  mean  plasma -density  gradients 
supports  the  concept  that  the  process  is  one  of  turbulent  cascade.  The 
damping  of  ESF  backscatter  at  low  altitudes  (below  250  km)  has  also  been 
noted  by  Tsunoda  et  al.  (1979)  and  Huba  and  Ossakow  (1981),  supporting 
the  damping  effects  of  electron-neutral  collisions  on  the  LHD  instabil¬ 
ity. 

The  most  comprehensive  measurements  of  small-scale  ESF  structure 
were  made  during  the  Plumex  I  rocket  flight  in  1979.  Simultaneous  mea¬ 
surements  were  of  irregularity  structure  that  spanned  five  orders  of 
magnitude.  In  situ  measurements  Inside  a  high-altitude  plasma  bubble 
associated  with  ESF  backscatter  revealed  that  LHD  waves  were  damped 
(Kelley  et  al..  1982).  On  this  basis,  Kelley  et  al.  (1982)  argued  that 
small-scale  structure,  those  smaller  than  the  ion  gyroradlus,  must 
result  from  wave-wave  coupling  and  not  linear  excitation  of  the  LHD 
instability.  Huba  and  Ossakow  (1979)  have  also  proposed  that  3-m  irreg¬ 
ularities,  those  responsible  for  Jlcamarca  radar  backscatter,  are  not 
linearly  excited  by  low-frequency,  collislonal  drift  waves  that  are 
heavily  damped  by  ion  viscosity.  Both  results  suggest  all  small-scale 
irregularities  below  the  ion  gyroradlus  are  produced  by  wave-wave  cou¬ 
pling.  The  role  of  the  LHD  instability  would  then  be  to  reduce  the 
amount  of  damping  in  the  wavelength  regime  below  about  1  m. 

Interpretation  of  Plumex  I  results  should  be  viewed  in  perspective. 
The  measurements  within  the  hlgh-altltude  plasma  bubble  were  made  during 
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the  decay  phase  of  Che  backscatter  plume.  Gradient  scale  lengths  on  the 
order  of  50  m,  as  observed  by  Costa  and  Kelley  (1978a,b)  and  Kelley  et 
al.  (1982)  are  probably  fairly  common  during  the  growth  phase  of  plasma 
bubbles  and  backscatter  plumes.  Using  the  threshold  curves  published  by 
Huba  and  Oasakow  (1981),  we  find  that  the  LHD  instability  is  operative 
in  regions  where  the  plasma  density  is  on  the  order  of  10^  el/cm*^. 
Tsunoda  et  al.  (1982)  have  shown  that  large  backscatter  plumes  are  asso¬ 
ciated  with  plasma  densities  of  that  magnitude.  On  this  basis,  we  con¬ 
clude  Chat  LHD  waves  are  linearly  excited  during  Che  growth  phase  of 
backscatter  plume  (and  plasma  bubble)  evolution. 

Once  the  irregularities  are  formed,  cross-field  diffusion  primarily 
controls  their  decay.  Classical  diffusion  would  result  in  the  persis¬ 
tence  of  irregularities  for  several  hours.  ALTAIR  measurements  have 
shown,  however.  Chat  the  decay  of  the  backscatter  plume  is  relatively 
rapid,  proceeding  at  a  rate  of  5  to  18  dB  per  10  min  (Tsunoda,  1980).  A 
good  example  is  the  decay  of  the  backscatter  plume  in  the  bottom  panel 
of  Figure  IS.  The  plume  on  the  left  side  of  Chat  panel  is  also  pre¬ 
sented  in  Figure  16,  mapped  about  30  min  later.  The  decay  is  evident 
near  the  head  of  the  plume.  The  decay  process  is  further  illustrated  in 
Figure  20.  The  segment  of  the  plume  located  in  Che  topside  F  layer  dis¬ 
sipates  in  about  30  min. 

Huba  and  Ossakow  (1981)  considered  diffusion  effects  on  ESF  irregu¬ 
larities.  They  concluded  Chat  anomalous  diffusion  must  be  occurring, 
probably  through  drift-wave  turbulence.  This  process  is  believed  to 
operate  via  the  universal  drift  instability  or  the  drlft-dlsslpaclve 
Instability,  at  spatial  wavelengths  less  than  1  km,  but  greater  than  the 
ion  gyroradlus.  As  these  low-frequency  drift  waves  smooth  out  the 
existing  gradients,  Che  growth  rate  of  the  LHD  instability  is  reduced. 
The  high-frequency  drift  waves  are  themselves  Incapable  of  smoothing 
gradients.  Because  diffusive  damping  is  proportional  to  k  ,  dissipation 
of  backscatter-producing  Irregularities  can  proceed  via  classical  (or 
anomalous)  diffusion.  Anomalous  diffusion,  however,  does  not  appear 
capable  of  dissipating  km-scale  structures.  That  is,  even  with  anoma¬ 
lous  diffusion,  these  km-scale  irregularities  would  take  eight  hours  to 
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IV  DISCUSSION  AND  CONCLUSIONS 


Late-tlne  (several  minutes  to  many  hours  after  burst),  hlgh- 
altltude  nuclear  effects  are  of  serious  concern  to  transionospheric  com¬ 
munication,  navigation,  and  surveillance  systems  because  of  the  severe, 
widespread  propagation  disturbances  that  ensue.  To  mitigate  these 
effects,  systems  designers  must  use  much  higher  frequencies  or  very 
elaborate  coding  or  equalization  schemes  or  both.  Both  alternatives  are 
very  costly  and  greatly  Increase  system  complexity.  Sound  decision  mak¬ 
ing  among  the  available  options  clearly  requires  accurate  phenomenologi¬ 
cal  predictions. 

It  has  been  convincingly  demonstrated  that  knowledge  of  the  spec¬ 
tral  characteristics  of  the  strlations  is  sufficient  for  predicting  the 
deleterious  effects  of  propagation  disturbances.  Thus,  it  can  be  argued 
that  the  main  objective  of  the  phenomenological  program  is  to  character¬ 
ize  the  temporal  and  spatial  variations  of  the  strlatlon  spectrum. 

Within  this  broad  objective,  however,  some  aspects  of  the  spectral  char¬ 
acterization  are  much  more  Important  than  others. 

The  spectral  strength  in  the  scale-size  regimes  that  dominate  the 
propagation  effects  is  most  Important.  For  frequencies  at  UHF  and 
below,  kilometer-scale  and  larger  structures  are  most  important.  At  GHz 
frequencies,  subkilometer-scale  structures  become  increasingly  more 
important.  It  should  be  kept  in  mind  that  under  strong-scatter  conditions, 
structures  larger  than  the  Fresnel  radius  dominate  the  scintillation 
effects.  Of  secondary,  but  not  necessarily  minor.  Importance  is  the 
spectral  shape  within  the  dominant  scale-size  regime.  Any  well-defined 
spectral  cutoff  is  clearly  important  insofar  as  extrapolation  to  higher 
frequencies  is  concerned. 


The  approach  that  ts  being  used  currently  in  nuclear  predictive 
codes  is  as  follows: 


(1)  Determine  where  in  a  HANE  environment  convective  instabilities 
are  strongly  driven. 

(2)  Estimate  the  maximum  integrated  structure  level  that  can  occur 
in  the  locally  unstable  volume. 


(3)  For  a  specified  outer  scale  wavenumber,  distribute  the  spec¬ 
tral  strength  according  to  a  k~^  power  law,  down  to  a  speci¬ 
fied  inner-scale  cutoff. 


This  procedure  establishes  a  reasonable  worst  case  based  on  our 
knowledge  gleaned  from  both  experiments  and  theory.  On  the  other  hand, 
the  approach  depends  critically  on  cutoff  wave  numbers  that  we  do  not  as 
yet  fully  understand.  The  ramifications  of  initial  conditions  are 
Important  as  are  the  diffusive  effects  that  control  the  evolution  and 
decay  of  the  structure.  Thus,  continued  research  is  an  Important 
adjunct  to  our  current  predictive  capability. 

As  discussed  in  Section  III,  measurements  of  ESF  phenomena  have 
provided  the  most  detailed  characterization  of  structure  development  in 
all  the  strlation  scale-size  regimes.  Indeed,  ESF  is  very  regular  in 
its  occurrence  and  characteristics;  moreover,  the  most  Intense  naturally 
occurring  irregularities  are  found  there.  Thus,  for  future  applica¬ 
tions,  ESF  phenomena  can  serve  as  a  systems  test  bed  to  demonstrate 
system  performance  in  a  well-defined,  reasonably  predictable  structure 
environment . 

The  primary  unresolved  issues,  such  as  the  effects  of  coupling  and 
large-scale  convection  patterns,  can  best  be  resolved  by  high-latitude 
measurements.  Moreover,  as  discussed  in  Section  II,  a  rich  variety  of 
other  potentially  relevant  phenomena  occurs  in  the  auroral  zone,  such  as 
energetic  particle  precipitation,  field-aligned  currents,  and  parallel 
electric  fields.  Finally,  current  phenomenology  models  cannot  make 
reliable  HANE  predictions  at  high  latitudes.  Thus,  future  effort  will 
emphasize  high-latitude  phenomenology  heavily. 

To  summarize,  current  HANE  phenomenological  codes  made  reasonable 
worst-case  predictions  of  propagation  effects.  Future  research  is 
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needed,  however,  to  refine  and  extend  the  predictive  capability,  particu¬ 
larly  where  very  expensive  mitigation  schemes  are  being  considered. 
Naturally  occurring  striations  provide  the  most  readily  accessible  data 
resource  for  verifying  our  ability  to  accurately  describe  the  relevant 
physical  processes  and  their  interaction  in  a  large-scale  system. 

The  research  reviewed  in  this  report  has  provided  both  an  extensive 
data  base  and  a  theoretical  framework,  within  which  the  relevant  physical 
processes  can  be  identified  and  extrapolated  to  refine  HANE  predictions. 
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